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Abstract 

We present an unifying description of the graded SL(p,q) KP-KdV hierarchies, in- 
cluding the Wronskian construction of their tau-functions as weh as the coefficients of 
the pertinent Lax operators, obtained via successive apphcations of special Darboux- 
Backlund transformations. The emerging Darboux-Backlund structure is identified as 
a constrained generalized Toda lattice system relevant for the two-matrix string model. 
It allows simple derivation of the n-soliton solutions of the unconstrained KP system. 
Also, the exact Wronskian solution for the two-matrix model partition function is 
found. 
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1. Introduction 



The Kadomtsev-Petviashvili (KP) hierarchy of integrable sohton nonhnear evolution 
equations 0] is among the most important physically relevant integrable systems. One 
of the main reasons for the interest towards the KP hierarchy in the last few years origi- 
nates from its deep connection with the statistical-mechanical models of random matrices 
((multi-)matrix models) providing non-perturbative discretized formulation of string theory 
0. Most of the studies in the latter area follow the ideas of the original papers 0], where 
the integrability structures arise only after taking the continuum double-scaling limit. There 
exists, however, an alternative efficient approach ^ for extracting continuum differential in- 
tegrable hierarchies from multi-matrix string models even before taking the continuum limit. 
More precisely, it is various reductions of the full KP hierarchy [constrained KP Hierarchies) 
which play the major role in the latter context. 

On the other hand, constrained KP hierarchies arise also naturally in purely solitonic 
context as shown below in subsection (6.2) (see also 0] and references therein). 

It is the aim of the present note to study various properties and provide exact solutions 
for a specific class of constrained KP hierarchies - the graded SL{p, q) KP-KdV hierarchies, 
which are intimately related with the two-matrix string model (which is the most physically 
relevant one). The following main results are contained in the sequel: 

(1) We establish the equivalence between conventional "symmetry" -constrained KP hierar- 
chies 1^ and multi-boson reductions of the full KP hierarchy p, |^, also known as graded 
SL{p, q')-type KP-KdV hierarchies 0, which appear in two-matrix models of string theory 
1^. In particular, we provide the explicit Miura map relating the above hierarchies. 

(2) Explicit exact solutions are found for SL{p, q) KP-KdV integrable systems, including 
eigenfunctions and r-functions, via special Darboux-Backlund(DB) transformations. 

(3) We establish equivalence between the set of successive DB transformations on the SL{p, 1) 
KP-KdV system and the equations of motion of a constrained generalized Toda lattice model, 
which embodies the integrability structure of two-matrix string models. 

(4) As a byproduct of (3) we obtain exact solution for this constrained Toda lattice system 
under specific initial conditions, relevant in the context of the two-matrix string model, and 
derive the exact expression for the partition function of the latter. 

(5) The present DB formalism provides a simple systematic way to obtain the n-soliton 
solutions for the full (unconstrained) KP system. 

2. Background on Generalized KP-KdV Hierarchies and Darboux- 
Backlund Transformations 

2.1 Constrained KP Hierarchies 

We shall consider the general class of constrained KP Lax operators with higher purely- 
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differential part also known as A^-generalized two-boson KP Lax operators |TD| Q : 

N r-2 N 

E 

j=i 



L = L+ + Y,<!>iD-'^i = D' + J2uiD' + Y,ai{D-k)-' 



(1) 
(2) 



One can also define an alternative consistent Poisson reduction of the standard KP hierarchy 
based on the pseudo-differential Lax operators of the form [§, 01 : 



N N N N 

L^=L^+Y. n n = ^+ + E n - <0 



-1 



i=l k=i i=\ k= 



(3) 

N-1 (4) 



which we will call multi-boson reduction of the full KP Lax operator. The above multi-boson 
reductions of the full KP Lax operators @ are defining the generalized graded SL{r + N, N) 
KP-KdV hierarchies pertinent to the string two-matrix models (cf. refs.[^, ^]). 

In [|l^] these two formulations of the constrained KP hierarchy have been related via 
successive DB similarity transformations. Below in section 3 we will establish their complete 
equivalence showing how the pseudo-differential Lax operators (|I]) and (0) can be rewritten 
into each other via a generalized Miura transformation. Due to this result we can limit 
ourselves to study DB transformations within the framework of the constrained KP hierarchy 
defined by A^-generalized two-boson KP Lax operators as in (|I]) 

2.2 On the DB transformations of the A^-generalized two-boson KP Lax operators 

The general form of a DB transformation on the A^-generalized two-boson KP Lax oper- 
ator (|1|) reads 0, : 



$0 
^0 



>i 



D 



-1 



X 



N 



(5) 
(6) 
(7) 



i=l 



X 



X 



N 1 

(x-'L+x) + E {x''^^) d-' i^^x) + <^'i*.) = {xDx-'l) x (8 



i=l 



$i = Xd. {x~'^i) 



x-'d-' {%x) 



where all functions involved are (adjoint) eigenfunctions of L (|l]), i.e., they satisfy: 



_d_ 

dtr, 



f = Llf f = X,<^^ 



_d_ 

dtn 



(9) 



(10) 



■^In order to avoid confusion, D will denote differential operator in the sense of pseudo-differential calculus, 
whereas derivative of a function will be denoted as d^f- 
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Let us particularly stress that the above eigenfunctions are not Baker- Akhiezer eigenfunctions 
of L from ([^) , unlike the construction in ref . |12 . 



We are interested in the special case when x coincides with one of the original eigen- 
functions of L, e.g. X = ^i- Then $i = and the DB transformation (|^) preserves 
the A^-generalized two-boson form (|l]) of the Lax operators involved, i.e., it becomes auto- 
Backlund transformation. Applying successive DB transformations in this case yields: 



^(fc+i) ^ . (,) ^ vl/f +1) = ($f , A; = 0,1,... (12) 



v^f+i) = -($W)-^9-(v^f)$f)) , . = 2,...,iV (14) 
Using the first identity (|TT]), i.e., L^^+^^T^^) = T^^^L^^^ , one can rewrite ([T^ in the form: 

^(fc) _^(fc-i)y(fe-2)...y(o) (^(^iio)y (15) 

whereas: 

^(fc) ^y(fc_i)y(fc-2)...y(o)$(o) ^ t = 2,...,N (16) 

Finally, for the coefficient of the next-to- leading differential term in (|I]) Ur-2 = r Res Lr = 
rS^lnr, we easily obtain from (^) (with x = ^i) its k-step DB-transformed expression: 



I {u^^ - 4%) = dl In ^ = 5^ In . . . (17) 



2.3 Wronskian Preliminaries 

Firstly we list three basic properties of Wronskian determinants. 

1) The derivative T>' of a determinant V of order n, whose entries are differentiable functions 
can be written as: 

D' = + D(2) + . . . + (is; 

where is obtained from D by differentiating the entries in the i-th row. 

2) Jacobi expansion theorem: 



W, if) Wu-i = WuWl_, if) - W[Wu^^ if) or d = ' (19) 



where the Wronskians are Wk = Wklipi, . . . jipk] and Wk-i (/) = . . . , V'fc-i; /]• For 

proof see Take a special class of Wronskians Wn = Wn[ip, . . . , d^~^ip] . Hence, from 
(H) we get : 



WX - WX = Wn W^^i (9» - C Wn^l (5» = Wn-lWn+1 ^ lu Wn 



(20) 
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3) Iterative composition of Wronskians: 



nn.,...T,{f) = (21) 

vVk 



where 



The proof of (|2T|) follows by simple iteration of ( |l9l) (see also the standard references on this 
subject [|15|, |16|, Q). For future use let us rewrite ( pTD as: 



(D + t;fc)(D + i;fc_i)...(Z} + ^;i) / = -^ ; t;. = 5,.ln-^ (23) 



2.4 DB Solutions of Two-Bose KP System and Connection with Ordinary Toda 
Lattice 

The two-boson KP system defined by the Lax operator L = D + $D^^\E' = D + a {D — b)~^ 
is the most basic constrained KP structure. We start with the initial "free" Lax operator 
L^^) = D and perform a DB transformation: 

LW = (^$(0)Z}$(0)"') D ($(o)Z}-i$(o)"') = D+ [$W (ln$W)"] D-' (24) 



The construction below is a special application of property 3) and eq. (^ll) . 

In the case under consideration the relevant formulas for successive DB transformations 
(IDjIO) specialize to: 

^i.k+i)^^{k) ^(^^(k)\j^^(k) ^ ^{k+i) ^ (^^(k)y^ ^26) 

whereas eq.(|l5|) acquires the form (proved easily by induction in k) : 
Introduce now 

0fc = In ^^^'^ = 6'^'= A; = 0, . . . (28) 

which allows us to rewrite ( PBD as a (ordinary one-dimensional) Toda lattice equation: 

5^0^ = 6'^'=+!-'^'= - e'^"-'*"'-^ (29) 

Introduce now new objects ipn as follows: 

0„ = ^„+i-Vn ^ ^„ = lnW^4$,a$,...,9"-i<I>] (30) 
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From equation ( ^0]) we find immediately an equation for ip^'- 



(31) 



with = for n < 0. We recognize at the right hand side of ( pTD a structure of the 
Cartan matrix for An- Leznov considered such an equation with Wronskian solution (in two 
dimensions) in |17]. 

Hence, the solutions of the (ordinary one-dimensional) Toda lattice equations, with 
boundary conditions ipn = for n < 0, reproduce the DB solutions of ordinary two-boson 
KP hierarchy (^) upon taking into account that $ = = exp (0o) = exp {^i). 

3. Equivalence between N- Generalized Two-Boson and 2N-Boson 
KP Hierarchies 

First, let us consider the simplest nontrivial case = 2 in ([1|). Applying the simple identity: 



\ A / 



w[xM ( w[x] j^.,w[xM' 



w[x] \w[xM' 



w[x] 



X-'D-\) (32) 



for arbitrary functions (j),ip,X) where in the second equality Wronskian identity (^l]) was 
used, we obtain: 



L+ + A^o^ (D-B. 



>{2) 



+ 4'^ [D-B, 



(2) 



D-B, 



(2) 



4^) 



(2) 



$1^1 + $2^'2 

W^[vi>2,Vi>i] 



5f 



(2) 



-dx In 



(33) 
(34) 

(35) 



Using successively the same type of identity as (p^), together with (pT]), e.g., for arbitrary 
functions 0, V'; • 

^ ^^V^ w[x] \w[x,uj^ i^fvi i 



+< 



W[x] _^-iW[x,uj] 



W[x] 



X-'D-\ 
X-'D-'X 



(36) 



we can prove by induction in that the A^-generalized two-boson KP Lax operator can be 
transformed into 2A^-boson KP Lax operator: 



r-2 



N 



N 



1=0 1=1 i=l 

= L, + j: A^r {D - )) {D - Bf^l) -'...{D-B 

i=l 



(N) 
N 



(37) 
(38) 
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upon the following change of variables, i.e., generalized Miura transformation: 

s=l ^ V^N, • • • , Wfc+iJ 

^(.) ^ -a.ln ^[^^^--^r'^-] (40) 

Let us now illustrate the equivalence between (|37D and (^) in the inverse direction. To 
this end it is more convenient to use the (r, q) form of (^. Let us define the quantity 

Qfc, = (-1)*-'= Jq,J /■■■J<lk {dx'y-'+' l<k<i<N (41) 

Then using that D^^Qi i^iqi = D^^Qi iD — we obtain from the first eq.(|^) 

N N 

Ln = L+ + Y. n D-\k + nD-^ i-Qi,N-iqN) (42) 

i=2 k=i 

where 

rf) =r, + rigi,,_i ^ = 2, . . . , iV (43) 
The above process can be continued to yield expression (|l|) with 

i—l i—k i—k i—k 

= fi ~\~ ^ "^k ^ ■ ■ ■ ^ ^ Qfc,i-Si_s,_i-lQi-s,;_fc_i,i-Si_fc_2-l 

fc = l Si_fc_l=Si_fe_2 + l S2 = Sl + lSl=l 

■ ■ ■ Qi-s2,i-s^-iQi-s^,i-i l<i<N (44) 

^Ar = gAr , = (-1)^" V / 9^-1 / • • • / 9^ (t^a:')^"* l<z<iV-l (45) 



4. Exact Solutions of SL(p,q) KP-KdV via DB Transformations 

Armed with the Wronskian identities from section 2.3, we can now represent the fc-step 
DB transformation ( [T5| ) — (|1^ in terms of Wronskian determinants involving the coefficient 
functions of the "initial" Lax operator 

r-2 N 

Lio) =D' + Y. ^?D' + ^ $f ^D-^M/f (46) 

1=0 1=1 

only. Indeed, using identity (pT]) and defining: 

A; = 1,2,... (47) 

we arrive at the following general result: 

Proposition. The k-step DB-transformed eigenfunctions and the tau-function ( ff^ j — ( [I^ j 
of the SL{r + A^, A^) KP-KdV system P% ) for arbitrary initial L^^^ are given by: 

= ""^ ''"'^ ^ , j = 2,...,Ar (49) 

Wk[^f\x^'\...,x^''-^)] 

r« = l^,[$S°),x^i),...,x('=-i)]rW (50) 
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where t^^\t^''^ are the r-functions of L^^\ L^^'^ , respectively, and x''*-* is given by (^^. 

As an example let us consider the SL{3, 1) KP-KdV Lax operator, i.e., r = 2, = 1 in 
(0) (the latter is pertinent to the simplest nontrivial string two- matrix model [|TT|) : 



L = D'^ + U + A{D - By^ = D'^ + u + <I>D-^^ 



(51) 

From the basic formulas for successive DB transformations (|l2D-(|lT]), applied to (|5lD , we 
have: 



lik) = 1)2 ^ ^(k) ^ ^(fc) ^-l^(fc) 



_^ l(k+l) ^ J-(k)i{k) frp(k) 



-1 



rp{k) ^ ^(fc) JJ f ^{k)\ ^ 



M^'^) = 2Res = 2dl In = 2dl In ■ ■ ■ <l>(°) 



9^ 



52^(fe-i) + 292ln f$('=-2)...$(0) 



<|)(A:-2) 



(52) 
(53) 
(54) 

(55) 
(56) 



«(o) = , = ; <l>(°) = / — c(A)e«(^'«) , ^ (A, {t}) = Ax + V \H, (57) 

27r fr^ 



i>2 



where in ( p7| ) is an arbitrary eigenfunction of the "free" L^^^ = (the contour T 
in the complex A-plane is chosen such that the generalized Laplace transform of c(A) is 
well-defined). 

As a corollary from the above proposition, we get in the case of (^2]) : 



$(fc) ^ . . . 2^(0) ^^2*:^ 



(0) 



W 



$(0)^52^(0)^ . . . 



W 



$(0),92$(0),...,9^('=-^)$(0) 



(58) 
(59) 



Substituting (|58|),(|59|) into 
coefficient functions of (FTD : 



we obtain the following explicit solutions for the 



= 292lniy [<l>(°),9,2$W,...,92("'i)<l>W 



w 


$(0),a2$(0),.. 


_ ^^2(n-l)^(0)" 


w 


<l>(0),a2$(0),.. 


.,a'("-2)$(0)- 



(60) 
(61) 

(62) 



Similarly, in the more general case of SL{r + 1, 1) KP-KdV Lax operator for arbitrary 
finite r : 



w 


$(°),9^$(°),...,(9^"$(o)' 


W 


$(0)^ ^2^(0)^ _ _ _ ^ ^2(n-2)^(0)" 


(w 







r-2 



«=0 



(63) 
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which defines the integrable hierarchy corresponding to the general string two-matrix model 
(cf. In], H), the generalizations of (|58|) and (^) read: 



(0) 



W 



w 



(64) 
(65) 



^ui% = Res l"- = dl In r^^^ , r^^) = [$(°), 9;<l>(°), . . . , Si'^-D-^W 
where is again given explicitly by (^). 

5. Relation to Constrained Generalized Toda Lattices 



Here we shall establish the equivalence between the set of successive DB transformations of 
the SL{r + 1, 1) KP-KdV system (||) : 



^(fc+i) ^ rj^ik) j^{k) (^{k)\^ 1 ^ rj^ik) ^ 



(66) 
(67) 



z=o 



and the equations of motion of a constrained generalized Toda lattice system, underlying 
the two-matrix string model, which contains, in particular, the two-dimensional Toda lattice 
equations. 

For simplicity we shall illustrate the above property on the simplest nontrivial case of 
^L(3, 1) KP-KdV hierarchy (|l|). We note that eqs.(||)-(|6D (or (|6b|)-(|2D) can be cast in 
the following recurrence form: 



(9^1nA("-i) 



with "initial" conditions (cf. (^)) : 

^(0) = ^(0) 



u 



(n-l) 



U 



W = 



= 9, In $ 



(68) 
(69) 

(70) 



(71) 



where $ is so far an arbitrary function. Now, we can view (|68| ) -([70| ) as a system of lattice 
equations for the dynamical variables A^'^\ B^'^\u'^'^'^ associated with each lattice site n and 
subject to the boundary conditions: 



U 



(n) 



, n<0 



(72) 



Taking (|7I]) as initial data, one can solve the lattice system (p8|)-(|70D step by step (for 



n = l,2. 



.) and the solution has precisely the form of 
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The lattice system (|0BD-([7nD can be identified with the ti = x evolution equations of the 
constrained generalized Toda lattice hierarchy defined as follows [HI 111] : 



d r r ■ 

g^Q = [Q[+),Q 



d - r ^ - 

—Q = 



1,...,P1 



-9 



Q,Q 



1 



(73) 

(74) 
(75) 



Here Q and Q are semi-infinite matrices, i.e., with indices running from to oo, with the 
following explicit parametrization: 



Q 



nn ^ flow 
Qnn 



Qn,n+1 — 1 ) Qn,n—k — 0,k{^) ^ 

for m — n > 2 , n — m > p2 



,P2 - 1 



Qnn = ho{n) 



Q 



n,n~l 



Q 



n,n+k 



h{n)R^li---R^lk A; = l,...,pi - 1 



Qr. 



for n ~ m > 2 



m 



n > pi 



(76) 



(77) 



The subscripts — /+ in ([73|),([7^) denote lower/upper triangular parts, whereas (+)/(—) de- 
note upper/lower triangular plus diagonal parts. In the case under consideration the number 
P2 = 3 in ([76|) , whereas the number pi in (^) is arbitrary finite or oo 0. 

Note the presence of the non-evolution constraint eg. (ffSD , which is called "string equa- 
tion". The lattice equations for the matrix elements ak{n) of Q (the first eqs.(^) and (|7i 
can be solved explicitly as functionals of the matrix elements of Q : 



P2-1 

Q{-) = E «-^(-) 

s=0 



is + l) 



t 



s+l 



(7J 



Furthermore, it is more convenient to introduce another matrix Q (with matrix elements 
Rn,bk{n), cf. (|77D) in place of Q defined as: 



Q 



P2-1 

(-) 



s=0 



(-) ' 



Q 



{+) 



Q 



P2-1 



Q 



(79) 



with I+nm — ^n+i,m, whcrc the last equality follows from (^). More generally, we have the 
relations Q'f^) = J2t=o IsaQI^) for any s = 1,. . . ,p2 with coefficients •jsa simply expressed 
through as (f^). Specifically we get: 

2 



Rr, 



7ii 



731 



lllRn 
(ttp2-ly 

3 

P2 - 1 



1 

P2-1 



bo{n) = 711^0 (^) + 
2 



721 



721 

2711 ' 

ttp2-2 



^i(^) = liih{n) + 



731 



37 



11 



721 



27 



^0) 



11 , 



-| po —3 



a 



P2-l> 



P2-3 



P2-4: 



a' 



P2-2 



P2-1 



P2-4 
I P2-1 



2(P2 



1 N 2p2-5 
^ («P2-l) ''^-^ 



^1) 



■^Both numbers pi.2 indicating the number of non-zero diagonals, outside the main one, of the matrices 
Q and Q are related with the polynomial orders of the corresponding string two-matrix model potentials, 
whereas the constant g in ( [75|) denotes the coupling parameter between the two random matrices. 
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Accordingly, the evolution eqs.([7iD acquire the form: 



-^Q = Q, Q'' with = 7s(jTr=^ 



12) 



The remaining independent lattice equations then read (we write down explicitly only 
the ti = X and ^2 evolution equations for the P2 = 3 case) : 



dxliaRn+i = bo{n + l) -bo{n) , bi{n) - bi{n - 1) = d^bo^n) 
Rn+i - Rn = (blin) + bi{n) +bi{n-l) 



d ^ r . . . 

TT^-Rn+l = dx dxRn+l + 26o(^)-Rn+l 
OT2 



(83) 
(84) 

(85) 



d_ 



kin) = dx 2h{n) + blin) - dxHn) 



d_ 



bi{n) = dxRn+i 



Now, we observe that the system of Darboux-Backlund equations for SL{3, 1) KP-KdV 
hierarchy (|6^)- ([70|) exactly coincides upon identification: 

= bo{n - 1) , = 26i(n - 1) , A^") = R^ (87) 

with the X = ti constrained Toda lattice evolution equations (p^)-(p^). Also, the higher 
Toda lattice evolution parameters can be identified with the following subset of evolution 
parameters of the SL{p2, 1) KP-KdV hierarchy ( p^ ) [|TT], ^ : 



. _ .KP-KdV , 
,P2 , Ir — lj.(^p^_'i^^ r—l,...,pi 



the second identification resulting from (ffOl). 

In particular, excluding bo{n) = and bi{n) = |'u("+^) in (|TDp using 

obtain the two-dimensional Toda lattice equation for A*^"^ = Rn '■ 



(88) 



, we 



dxd{^ In = - 2A(") + 



(89) 



6. Discussion and Outlook 



6.1 Partition Function of the Two-Matrix String Model 

The partition function Zj^j of the two-matrix string model is simply expressed in terms 
of the Q matrix element bi{N — 1) at the Toda lattice site — 1, where A^ indicates the size 
(A^ X A^) of the pertinent random matrices: O^IyiZn = bi{N — 1) (cf. 0, |ll|]). Thus, using 
(|65|) and (|57| ) together with (^), and accounting for the relations (pO|)-(|8TD, we obtain the 
following exact solution at finite N for the two-matrix model partition function: 



W 



$(0) 




exp — A^ 



/ |^c(A)e<^'^*'*» , i (A, {i, t}) ^ E A% + E 



(90) 



(91) 
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where x = ii and the 7-coefficients are defined in (KTI). The "density" function c(A) in (^Tj) 
is determined from matching the expression for ^^^'\dx In <l>^°) = 6o(0) = 7ii&o(0) + 721/2711 
(cf. ([7I|),(|87|),([80|)), with the expression for 60 (0) in the orthogonal-polynomial formalism [Q: 



^ exp |'y" ^ ^ c/AirfAa exp X[tr + ^ Aat, + fi-AiAa 



(92) 

tl=tp2-l{t) 



Obviously, the most important question now is to study the physical double-scaling limit 
0] of (pop , which amounts to a special fine-tuned limit 00. The latter presumably 

includes renormalizations and critical point approaching of the (tr,is^ parameters. 

6.2 Connection to Grassmannian manifolds and n-Soliton Solution for the KP 
Hierarchy. 

Let {ipi, . . . ,^|Jn} be a basis of solutions of the n-th order equation Lip = 0, where L = 
{D + Vn) (-D + Vn-i) . . . (-D + fi). If Wk denotes the Wronskian determinant of {ipi, . . . , ipk} 
then one can show that |18|, |T6[ : 

v, = d(ln^^ Wo = l (93) 

This allows to show that the space of differential operators is parametrized by the Grass- 
mannian manifold (see e.g. [0, [T^). Start namely with the given differential operator 
Ln = Z)" + MiD"~^ + - ■ ■ + Un and determine the kernel of L„ given by n-dimensional subspace 
of some Hilbert space of functions H, spanned, let say, by {ipi, ■ ■ ■ , 'ipn}- This establishes the 
connection one way. On the other hand let {ipi, . . . , V'n} be a basis of one point fl of Gr*-"-* 
being a Grassmannian manifold. Define the differential equation as Ln{fl)f = Wk{f)/Wk. 
From (|23| ) this associates the differential operator 

Ln = + UiD''~^ + ---+Un = {D + Vn){D + Vn-l) ■ ■ ■ {D + Vi) (94) 

given by a Miura correspondence to a given point of the Grassmannian. 

Recall now correspondence (equivalence) between iV-generalized two-boson KP and 2N- 
boson KP systems — eqs.(p7|)-(^). The relation (|40| ) has the form like in (|23|) and, therefore: 



D-BrY{D-Bi^}y'...{D-B^,'Y 



= ; i<j<N (95) 



The above relations generalize the relations encountered in the study of flags manifolds and 
clearly deserve more investigations. 

Let us now comment on connection to n-soliton solution for the KP hierarchy. Assume 
that the above functions ipi i = 1, . . . ,n have the property dmipi = d'^ipi for arbitrary m > 1 
(dm = gf-)) in other words tpi are eigenfunctions of L^") = D. We introduce L = LnDL^^, 
where L„ is deflned in terms of {V'l, • • • jV^n} as in (|93D and (^4]). It is known that such a 



Lax operator satisfles a generalized Lax equation dmL = [L™, L] [||, [22 . 
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Using ( pri) we can rewrite the above Lax operator as a result of successive DB transfor- 
mations applied on D: 

L = LnDL-' = Tn Tn-1 " " " Ti DTr' ■ ■ ■ T-_\ T'' (96) 



where Tj are given in terms of Wronskians as in (p2D. It follows that L can be cast in a 
form of the Lax operator belonging to the n-generalized two-boson KP hierarchy having 
the form as in (|^) with r = 1, N = n. Using the formalism developed in this paper one 
can prove by induction that the corresponding r-function of L takes a Wronskian form 
Tn = Wn[ipi, ■ ■ ■ yipn] reproducing n-soliton solution to the KP equation derived in In 
fact, choosing ipi = exp (j^tkO^i^ + exp (j^tkf^i^ allows to rewrite r„ in the conventional 
form of the ra-soliton solution to the KP equation 
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Note Added. After completion of this paper we became aware of refs.|23| where Wronskian 
expressions for partition functions of matrix models have been obtained by a different method. 
We would like to stress that our result (|90[)-(^),(p2|) explicitly incorporates the "string-equation" 
constraint (^) on the Toda lattice system. 
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